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AGGREGATION RATES IN ONE-DIMENSIONAL STOCHASTIC 
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St. Petersburg State University, Universite Lille I and Universite Paris VI 

We consider one-dimensional systems of self-gravitating sticky 
particles with random initial data and describe the process of aggre- 
gation in terms of the largest cluster size Ln at any fixed time prior 
to the critical time. The asymptotic behavior of L„ is also analyzed 
for sequences of times tending to the critical time. A phenomenon of 
phase transition shows up, namely, for small initial particle speeds 
("cold" gas) Ln has logarithmic order of growth while higher speeds 
("warm" gas) yield polynomial rates for Ln- 

1. Introduction. We consider a system of n particles living in one-dimen- 
sional space. At initial moment, every particle is characterized by its mass, 
initial position and initial speed. There exists a pairwise gravitation between 
the particles. Between the moments of shocks, particles move in this grav- 
itation field according to usual rules of Newtonian mechanics. The shock 
between two (or more) particles results in the birth of a new particle ( "clus- 
ter") whose characteristics are defined by the laws of conservation of mass 
and momentum (while the energy is dissipated at these nonelastic collisions) . 

The clusters gradually become larger and larger while the number of clus- 
ters diminishes — until the unique cluster, containing the totality of mass, 
remains on the line. This collapse vaguely models emergence of a "star" 
from dispersed "dust." Indeed, the roots of the model are in astrophysics 
(see [14, 18]), but the reason of actual interest in similar particle systems is 
due to their relation with solutions of nonlinear PDEs such as the Burgers 
equation (see [1, 3, 4, 7] and the references therein). 

The aim of the present work is to describe essential features of the aggre- 
gation process provided that initial data are random and that the number 
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of initial particles tends to infinity. Assuming that at the beginning there 
exist n particles, we try to understand the behavior of Ln{t) — the largest 
number of particles that form a common cluster at time t. Therefore, Ln{t) 
increases, as time goes by, from 1 at the beginning to n at the collapse time. 

We will extend the quantitative results obtained in recent works [[2, 7, 10], 
[16, 17]]. It is known from these works that, under reasonable assumptions, 
there exists a critical time T* such that lim„_>oo = in probability for 
every t <T* and lim^^oo -^^^i^l = i probability for every t>T*. In other 
words, the time when essential collapse occurs, is, in fact, deterministic. 
"Essential" means that after critical time one can still observe a number of 
small peripheral clusters, but the main core has already been formed and it 
contains the overwhelming part of the total mass. 

Our results show that the aggregation process behaves rather differently 
in the cases of "large" and "small" initial speeds. In the literature the most 
frequently used model is that of i.i.d. initial speeds. For this case, which is 
naturally interpreted as "warm" system (see, e.g., [10, 17]), we essentially 
show that for any t <T* , 

L„(t)«c(t)n2/3(logn)V3. 

On the other hand, whenever initial speed is small or vanishes, which is 
natural to interpret as "cold" system (see [7]), one has 

Ln{t) ^ c{t)logn. 

Our results, in fact, cover the whole range of possible behaviors containing 
the two aforementioned special cases. 

In any case, one can observe from these formulae that L„(t) <C n, that is, 
the aggregation process is rather slow and that genuinely macroscopic clus- 
ters appear only shortly before the critical time. Therefore, it is interesting 
to describe the behavior of L„(t„) for t„ — > T*. This is done for cold systems 
in Theorem 3.4. 

We give a rigorous description of the model in Section 2 and state our 
results in Section 3. Section 4 provides the necessary information on con- 
servation laws which control the behavior of the systems. Finally, Section 
5 contains the collection of proofs. There is an amazing contrast between 
the elementary nature of the results and the rather advanced techniques one 
needs to obtain them. 

2. Systems of sticky particles. 

2.1. Dynamics of deterministic systems. We consider the following one- 
dimensional system (gas) of particles with gravitation. At the starting mo- 
ment t = 0, our gas consists of n particles, positioned on the real line at 
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points with coordinates xi{0), . . . ,x„(0). Particles are always enumerated so 
that xi(0) < • • • < Xn{0). Each particle is characterized by its mass rrii and 
its initial speed Vi{0). The particles move under the action of forces of pair- 
wise mutual gravitation. The gravitation force acting on a particle of mass 
m positioned at x, from a particle of mass fi positioned at y, is 

(2.1) F = jmij,sign{y — x). 

Here, 7 > is a positive gravitation constant. Note that there is no depen- 
dence of the gravitation force on the distance between the particles. This 
feature is typical for one-dimensional gravitation models. 

Between collisions particles obey the Newton's second law {F = m^^). 
The total gravitation force acting on a particle is obviously proportional to 
the difference of the masses on its right and on its left. 

On collision particles stick together, following the conservation laws of 
mass and momentum. In other words, two particles with characteristics 
(m, v) and (//, w) produce one particle with mass M = m + n and speed 
V = {mv + ^w)/M. This is a completely nonelastic collision, kinetic energy 
is dissipated. The particles which are born upon collisions are called clusters. 

We assume that the initial particles do not die at collisions but continue 
their movement as parts of the created clusters. Therefore, the position, 
the speed and the acceleration of a particle are understood as those of the 
cluster containing this particle. We denote them by Vi{t) ■.= x[{t) and 

x'-{t), respectively, for each t>0. The destiny of each particle is therefore 
defined during the life time of the system. The speed and the acceleration 
have jumps at the times of collisions of the particle or of the cluster which 
contains it (actually, the acceleration also has 5-s at collision times). 

Throughout the paper, we call particles only the initial particles, and 
clusters the products of collisions as well as the initial particles. Thus, at 
any time, the system consists of a number of clusters, each cluster being a 
set of one or more particles. 

When Vi{0) = (zero initial speed), the system is referred to as a cold gas. 
If a gas is not cold, we call it a warm gas. 

2.2. Similarity of systems. Let positive numbers (similarity coefficients) 
Cx, Cv, Cm, ct and Cy satisfy the conditions 

^ Cx Cy 

CyCt CyyiCyCl 

Consider two systems of particles. The first system has initial data (xi(0), Ui(0), mj) 
and gravitation coefficient 7, while the second has initial data {cxXi{fi),CyVi{Q),Cm'mi) 
and gravitation coefficient 0^7. These two systems are similar at any time 
t > 0, that is, the first system at time t has at position x a cluster of mass 
m moving with the speed v if and only if the second system at time ctt has 
at position CxX a cluster of mass Cmm moving with the speed CyV. 
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2.3. Stochastic systems of particles. Let us introduce a stochastic version 
of the problem by considering initial parameters of the particles as random 
variables, while the interaction rules and the dynamics remain deterministic 
and follow the rules described above. Moreover, we consider the asymptotic 
situation, in which the number n of particles tends to infinity, while other 
parameters vary in a reasonably consistent way. 

2.3.1. Masses. We assume that the masses of the particles are determin- 
istic and equal to pn~^ each, where ;9 > is a fixed constant. This leads, at 
the limit, to a model with a given density of the matter p. We call the size 
of a cluster the number of particles in it. 

2.3.2. Initial speeds. We assume that the initial speeds of the particles 
have the form Vi{0) = anUi,l <i<n, where (uj) is a collection of i.i.d. 
random variables with zero mean and unit variance. Zero mean is assumed 
only for convenience of notation. It is easy to see that adding the same 
constant to all initial speeds affects neither the times of collisions nor the 
sizes of the created clusters. The only consequence is an additional uniform 
drift of the whole system. 

In the literature, the scaling parameter an is usually assumed to be in- 
dependent of n. Then the case (T„ = corresponds to a cold gas, while 
(7„ = o" > corresponds to a warm gas. We prefer to handle the setting with 
variable an, since we are interested in identifying the border separating the 
models whose properties are close to those of cold and warm gas, respec- 
tively. 

2.3.3. Initial positions. Three meaningful different models are distin- 
guished here, though they often lead to the same asymptotic results. 

The lattice deterministic model assumes that the particles are initially 
located on the lattice Xi{0) = i/n, 1 <i <n. 

The Poisson model assumes that the particles are initially located on the 
positive half-line at the first n points of a Poisson point process with in- 
tensity n. In other words, they are located at the times of first jumps of a 
Poisson process of the just mentioned intensity (the space where the par- 
ticles live is interpreted here as a time parameter of the process). By the 
well-known property of Poisson processes, the differences Xi[{)) — Xj_i(0) [no- 
tation: xo(0) := 0], 1 < i < n, are independent exponential random variables 
with mean n~^. 

The i.i.d. -model assumes that the particles are initially located at the 
points corresponding to a sample of n independent random variables uni- 
formly distributed on the interval [0, 1]. Since the particles are indexed by the 
order of initial positions, the initial location of the ith particle corresponds 
to the ith order statistics of the uniform sample. 
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2.3.4. Relation between the Poisson and the i.i.d. -models. The Poisson 
model is more convenient for investigation due to the (aforementioned) prop- 
erty of independence of distances between the particles at time zero. The 
following passage from the Poisson model to the i.i.d. -model is useful; it is 
well known and stated here without proof. 

Fact 2.1. Fix fc > 1, and let Xi, . . . , X^, X^+i be the times of first k 
jumps of a Poisson process of constant intensity. Then the random variables 
Xi = 1 < "i < ^, have the same joint distribution as the order statistics 

of a sample of k i.i.d. random variables uniformly distributed on the interval 
[0,1]. Moreover, the random vector and the random variable 

Xj^j^i are independent. 

Note that the intensity of the Poisson process does not play any role. In 
the sequel we use a standard Poisson process (of unit intensity), as well as 
a Poisson process of intensity n. 

3. Main results on aggregation rates. In this section we investigate the 
asymptotic behavior of Ln(t), the largest cluster size at time t in the system 
of n initial particles. Since the cluster size coincides with the number of 
initial particles in it, we always have 1 < Ln{t) < n. 

The collapse time T[n) is the first time when all the particles belong to 
a common cluster, that is, when the total collapse of the particle system 
occurs. 

Let us introduce the critical time of the system defined by 
(3.1) T*:={^p)-^l\ 

3.1. Aggregations in a cold gas. First we recall a known result (see, 
e.g., [7]) on the collapse time, which motivates the definition (3.1). 

Fact 3.1. In a cold gas, for any of the three models (lattice, Poisson or 
i.i.d.) of initial positions, the random variable T{n) converges in probability 
(when oo) to the deterministic constant T* . 

Therefore, the behavior of the largest cluster size after the critical time is 
trivial: 

Corollary 3.2. In a cold gas, for any of the three models (lattice, 
Poisson or i.i.d.), for any t > T* , we have lim„_^oo IP{-^n(i) =n} = 1. 

We consider now L„(t) for t <T* , that is, we are interested in the largest 
cluster size at a fixed time prior to the critical time. 
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The case of lattice initial positions is particularly instructive. For a cold 
gas we deal here with a purely deterministic model. All the particles move 
to the point of general meeting without any collision. More precisely, the 
particle i moves along the trajectory 

x.{s) = - + ^{n-2i + l), 
n 2n 

with constant acceleration. We see that at time T* all the particles simul- 
taneously meet at the barycenter point x = . The absence of collisions 
prior to the global meeting follows, for example, from the formula 

Xi{s) - Xi-l{s) = . 

n 

However, in other models of initial positions, the situation is not so trivial. 

Theorem 3.3. In a cold gas for both Poisson and i.i.d. -models, for any 
t <T* , we have 

(3.2) lim ^^^ = /((t/r*)2)-^ in probability, 

n^oo logn 

where 

(3.3) /(r) := r — 1 — logr, r > 0. 

The logarithmic order, already established in [7], indicates a very slow 
growth of clusters in a cold gas. Our contribution is to prove the existence 
of the limit and to determine its exact value. 

The function /(•) in (3.3) is the rate function in the large deviations for 
the exponential law. 

We set aside the delicate question of studying the behavior of L„(T*), 
and consider the maximal cluster size in a cold gas at times shortly prior to 
the critical instant T* . 

Theorem 3.4. In a cold gas for both Poisson and i.i.d. -models, for any 
sequence (t.„) C [0,T*] which satisfies assumptions 

(3.4) lim tn = T* and lim n(r* — 1„)^ = oo, 

n— >oo n— >oo 

we have 

y Ln{tn){T*-tn? {T*f . 

hm — -Tf- = m probability. 

n^oo log(n(r* - tnY) 2 F y 

Assumption (3.4) is quite natural, since it is shown in [7] for total collapse 
time T(n) that T* — T{n) has order n"^/^. The assertion of our theorem ob- 
viously becomes false when (t„) approaches T* faster than what is admitted 
in (3.4). 
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3.2. Aggregations in a warm gas. Consider now a warm gas, that is, a 
system of particles with nonzero initial speeds Vi{0) = (T„nj, where the speed 
scale ((T„) is a sequence of nonnegative numbers and (uj) is a family of i.i.d. 
random variables with E(uj) = and lE(n^) = 1. 

The aggregation in a warm gas has been extensively studied in [2, 10, 
16, 17]. The critical time T* is still defined by (3.1), but its interpretation 
becomes slightly different. It is not interpreted any more as a limit value 
of collapse times, but as a time after which there exists a single cluster 
of huge mass and, possibly, a dust of small clusters whose total mass is 
negligible. In other words, analogously to Corollary 3.2, for any t > T* , we 
have lim„_>oo ^"^^^ = 1 in probability. 

The aggregation in a warm gas at critical time T* is remarkably studied 
in [17] (for lattice model with (t„ = 1). Therefore, it remains to evaluate 
L„(t) at times t < T* (cf. Theorem 3.2.1 in [17]). 

Before stating the result, let us explain a natural but somewhat nonstan- 
dard notation of superior and inferior limits in probability which will be 
frequently used in the sequel. Let be a sequence of random variables and 
c G M. Then 

lim inf > c in probability 

n— ♦oo 

means that for all e > 0, 

hm P(^„ < c - e) = 0. 

n— >oo 

Similarly, 

lim sup < c in probability 

n— >oo 

means that for all e > 0, 

lim P(Cn>c + e) = 0. 



Theorem 3.5. Assume that 

(3.5) lim — ^ = oo and lim — ^ = 0, 

n^oon^^logn n^oo 72^/^ 

and that Ui has a finite exponential moment Eexp{a|iij|} < oo for some 
a > 0. Then, for any of the three models (lattice, Poisson or i.i.d.), for any 
t <T* , we have 

(3.6) liminf- No/sn"'^!'' / 2\Mrt - '^i(^) in probability, 

n^oo [nan) '^[log{n/a:^))^'-^ 

(3.7) limsup ^^^-^ — < C2(t) in probability, 

(na„)2/3(log(W'))V3 
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where 

/ f \ 2/3 

(3.8) c,{t):=2[^-^——^j , C2(t):=(20/3)i/3ci(t). 

In the basic case an = ^, we get the estimate Ln{t) ~ n^/'^ (logn)^^^. 

For relatively large an (high temperature gas) , a critical order turns out to 
be an ~ When approaching this order, the main term in the asymptotic 
expression for Ln{t) has order n, while the logarithmic terms of our upper 
and lower bounds are not of the same order anymore. This would indicate 
a kind of phase transition. The critical time for such a hot gas should not 
follow the formula T* = (7p)~^/2, but tend to infinity. 

It is interesting to understand what happens if we replace the assumption 
Eexp{a|iii|} < oo by the weaker condition 

(3.9) E[\ui\P]<oo. 

Then a small number of particles with high initial speeds can substantially 
perturb the behavior of Ln{t). However, the following result says that if p 
is large enough, then Theorem 3.5 still holds true, but in narrower zones of 
speed ranges. 

Proposition 3.6. Assume that (3.5) and (3.9) hold. Let£>0. 

(a) Whenever p> 2 and an ^ £n^'^~^^^^ , we have (3.6). 

(b) Whenever p> 6 and an> £n^'^~^'^^^^'~^\ we have (3.7). 

It is worthwhile to compare Proposition 3.6 with estimates in [17]. In our 
notation. Theorem 3.2.1 in [17] yields Ln{t) < V/i > for an = ^ 

and p > 4:, while Proposition 3.6 provides a better bound L„(t) < const x 
n2/^(logn)^/^. These bounds become closer when p— > oo. 

On the opposite side of the scale, for small cr„ (we call it a low temperature 
gas), namely, for (T„ ^ logn, it is easy to establish the same behavior as 
for the cold gas, that is, (3.2) is true. 

Moreover, if we formally apply Theorem 3.5 to critically small an = 
const X n~^logn, we also get Ln{t) ~ logn just as in Theorem 3.3. More 
rigorously, we have the following analogue of Theorem 3.5. 

Theorem 3.7. Assume that 

(3.10) lim — ^ = c G (0, cx)), 

n-*oo logn 

and that Eexp{a|iij|} < oo for some a > 0. Then, for any of the three models 
(lattice, Poisson or i.i.d.), for each t<T*, 

liminf /^^^ > Ci(t) in probability, 
n— >oo logn 
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limsup— — < C2(t) in probability, 
n-»oo log n 

where < ci{t) < C2{t) < oo are constants. 

4. Barycenter technique. In this section we describe and extend a barycen- 
ter technique for the study of particle systems. This technique has been 
already used, for example, in [2, 6, 7, 10, 17]. 

We identify the particles with their numbers, and call a block of particles 
any set of particles which have consecutive numbers. We denote such a block 
by J = {i,i + k]; this block contains the particles numbered from i + 1 to 
i + k. 

The number k is called the size of the block. It is worthwhile to mention 
that at any time t>0, some particles in the block may belong to a cluster 
containing particles that are not in the block. 

We say that a block is free from the right up to time t if none of its 
particles has collided up to time t with any of the particles initially located 
to the right of the block. A block free from the left is defined similarly. 
Finally, a block is free up to time t if it is free both from the left and from 
the right. We note that collisions inside a free block are possible. 

Let Mj := J2j£j^j- Define the barycenter of the block J by xj{t) := 
Mj^Y.j^jm,x,{t). 

Define 

(4.1) x*j{s) := xj{0) + x'j{0)s + lj{Mf - Mf^)s\ 

where Mj^"* and Mj^^ denote the total masses of the particles to the right 
and to the left of J, respectively. Note that s i— > x*j[s) represents the tra- 
jectory of the barycenter of the block J without taking into account the 
collisions with external particles and is completely expressed in terms of 
initial data of the particles in the block J. 

The following observation contains a basic idea in the study of such par- 
ticle systems: 

"The barycenter of a free block moves with constant acceleration, as if it 
were a single particle with mass equal to the total mass of the particles in 
the block." 

More precisely, we have the following. 

Proposition 4.1. Let a block J be free from the right up to time t. 
Then 

(4.2) xj{s)>x*j{s), se[0,t]. 
Similarly, if a block J is free from the left up to time t, then 



(4.3) 



xj{s)<x*j{s), se[0,t]. 
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The proof of Proposition 4.1 is straightforward, and is omitted. An im- 
mediate consequence is the following 

Corollary 4.2. Let a block J be free up to time t. Then 
xj{s)=Xj{s), 0<s<t. 

We call the collapse time t j of a block the first time when all the particles 
of the block J are in a common cluster. It is possible to express ij via 
the random variables of type x\;A C {1, 2, . . . , n}. Toward this aim, let us 
introduce for every particle j, the first time fj^ when it collides with the 
neighboring particle j + 1. It is clear that for J = {i,i + k], we have 

(4.4) t}i=maxttr- 

0<r<fe 

On the other hand, let us consider for each triplet a < j < f3, the quadratic 
function 



We know that 



Moreover, the main coefficient of the quadratic function is positive since by 
(4.1) it is equal to 

Therefore, there exists a unique j ^ > such that 

Qa,j,l3{T*jA = 0- 

From the mechanical point of view, r* ^ ^ is the time when the trajectories of 
barycenters of blocks {a,j] and (j,/3] meet, as long as we do not take shocks 
into account. The value T*jp has an explicit expression via the coefficients 
of the polynomial Qij.js, in terms of initial positions, speeds and masses of 
the particles. It is worth noticing that for any t > 0, 

(4.5) *><j,/3 ^ Qaj,/3(t)>0. 

It turns out that the following is true. 



Proposition 4.3. For any j < n, we have 

min - 



(4.6) ^ = .J!^A^,P- 
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Proof. At any time s <tf, the particles j and j ' + 1 belong to different 
clusters, so that any block of the form (a, j] is free from the right and any 
block of the form (j,/3] is free from the left. By Proposition 4.1, 

Therefore, Qa,j,f3{s) < 0. It follows that tf < oo and by continuity of Q, we 
have Qaj,(}{tf) < 0. Hence, r* ^- ^ > fj. Since a < j and (3> j are arbitrary, 
we have tf < min„<j,/3>j T*j^p. 

To prove the inequality in the other direction, we note that by definition. 



a <j and P> j such that the blocks (a, j] and (j, /3] are free up to time tf 



tj is the collision time between particles j and j + 1, so that there exist 

< j and (3 > j sucl 
(and collide at time tf). By Corollary 4.2, 



Hence, tf = t* ^. It follows that tf > minQ,<j^ f^yj T^^jp- The proposition is 
proved. □ 

When estimating the collapse time of a block, it is sometimes more con- 
venient to deal with simpler expressions. For example, it is a consequence of 
(4.4) and (4.6) that for J := + k], 

(4.7) t}' < max r* i , ^ g Va < V/3 > i + fe. 

0<r<fc ' 

In particular, taking a = i and (3 = i + k gives 

(4.8) tj < max T^+,,,+fc. 

0<r<fc 

Note that the expression on the right-hand side depends only on the initial 
data of the particles in the block itself. 
There is a situation when (4.8) is sharp. 

Proposition 4.4. // a block J = + k] is free up to its collapse time 
tj, then 

(4.9) tj = max T*i^r,i+k- 

0<r<k 

Proof. It is sufficient to prove the inequality opposite to (4.8). Take 
arbitrary integer r £ [l,k) and consider the blocks {i, i + r] and (i + r,i + k]. 
We know that, up to time tj, the first block is free from the left, while the 
second is free from the right. Therefore, by Proposition 4.1, 

^(j,j+r](*j) — ^{i,i+r]{tj) = X(i+r,j+fc] (* J ) > X(j_,_^^j_,_^] (t j ) . 

Hence, Qi^i+r,i+k{tj) > 0, so that by (4.5), tj > r*^^^^.,.^. Since r G [1, A;) is 
arbitrary, this yields tj > maxo<r<fc T*i^^^i^f^. □ 
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5. Proofs. 



5.1. Proof of Theorem 3.3. In Theorem 3.3, we work with the cold gas, 
that is, there is no initial speed. Thus, 



(5.1) 



2n 



1/2 



We will first provide a lower bound for Ln{t) for the Poisson model, that is, 
if the initial positions (xj(0), 1 < i < n) of the particles are the first n jump 
times, denoted by Xi < • • • < of a Poisson process of intensity n. Recall 
that 



(5.2) 



Lr, (t) = max< /c : min t'^}<t\. 



Let /c < n. Consider the blocks Jr.= {lk,{l + l)k], Q < l < ^ - I. It follows 
from (5.2) that 



(5.3) 



P(L„(t) <k)< P(mmt^^ > t 



For a generic block J = (i, i + A:], we have, by (5.1), 

K 

(5.4) 



h{J) := max T*i+ 

0<r<fe 



2n \ 



max 

0<r<fc 



1 



i+k 



i+r 



k — r.'^^ , 1 

J=«+r+l 



Taking J = Jg 



j=i+i 

z, {i + 1)A;], and in view of (4.8), we have 



1/2 



tj<b{Ji) 



2n\i/2 



max 

0<r<fc 



(5.5) 



1 



ii+l)k 



j=^fc+r+l 



ik+r 



j=^fc+i 



1/2 



with the notation Xq := 0. Note that the random variables h{J^)^ for < 
^ < 1^ — 1, are independent and identically distributed. Therefore, in light of 
(5.3), we obtain 



(5.6) 



¥{Ln{t) < k) <F(mmb{Ji)>t 



[1 - P(6(Jo) < t))" < exp{-z.P(6(Jo) < t)}, 
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where v := [^J denotes the number of blocks. According to (5.4), 



P(5(Jo) < i) =P max 



0<r<fc \ k — r . 

\ j=r+l 



2n 



At this stage, it is useful to recall from Fact 2.1 that 



(5.7) 



Since (Xi,...,Xfc) and X^+i are independent, we arrive at the following 
estimate: for any (5 > 0, 



P(6( Jo) < t) > P ( max ( ^— 
\0<r<fc\ k-r 



1 



1 



(5.8) 



j=r+l 



xPfXfc+i< 



(1 + 2S)n 

Consider the empirical quantile function 

(5.9) Qfc(s):=inf{t:Ffc(t)>s}, < s < 1, 

where -Ffc is the empirical distribution function based on the random vari- 
ables (Xj)i<j<fc, that is, Ffc(t) := : 1 < j < /c, < t\. Then 



= Qa:(s), 



and we have 

k 

i 

A; 



E ^,-jE^, = Y^/Qfc(^)rf«-^/Qfc(^)'i^, 



j=r+l 

Accordingly. 



r 



max 

0<r<fc 



(5.10) 



k — r 



" ~ 1 ~ 1 

E ^.-^E^.-^ 



j=r+l 



< max 
0<t<l 



f\Qk{s)-s)ds-l f\Qkis)-s)ds 

l—tJt t Jo 



The Glivenko-Cantelli theorem for quantile processes (see, e.g., [15], page 95) 
asserts that limfc^oosupg \ Qk{s) — s| = almost surely. Hence, 



lim ] 



>( maxf-^ V Xj--j2X,]<l + 5] 



1. 
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For the second probability expression on the right-hand side of (5.8), we 
note that (i^23)n ^^ ~ {i+2S}n i't/T*)'^. Therefore, for ah sufficiently large k 
(how large depending on 6), 

(5.11) 

>p(^!^±l < {l + 36y\t/T*fy 

By Chernoff 's large deviation principle, for all large k, 

lp(^^PY<(l + 3<5)-i(t/r*)2^ >exp{-(l + 5)/((l + 35)-i(Vr*)2)(fc + l)}, 

where the large deviation rate function /(•) of the exponential law is as in 
(3.3). We choose now 

I. i.( \ 



(l+3(^)/((l+3<5)-i(t/T*) 

so that by (5.8), 



i/P(6(Jo) < t) > const x-^n-(i+^)/(i+2'5) ^ oo. 

logn 

We derive via (5.6) that P(L„(t) <k)^Q. Therefore, 

liniinf > (l + 35)^i/((l + 3(5)^1 (Vr*)2)-^ in probability. 

n^oo logn 

By sending 6 to zero, we obtain the desired lower bound in Theorem 3.3 for 
the Poisson model: 

liminf > I{{t/T*fy^ in probability. 

n— »oo log n 

We prove now the upper bound for the Poisson model. Let k <n. Assume 
that we are in the situation Ln{t) > k, which means that prior to (or at) 
time t, a cluster (say, ^4) of size A; + 1 or larger appears. Actually, one can 
choose A of size less than or equal to 2k, but we do not need the upper 
bound in this proof. Consider the block of particles J which corresponds to 
^ in a natural way. Notice that J is free up to the time of formation of A. 
By (5.4) and (4.9), we obtain 

(5.12) b{J)=f}<t. 

To bound 5( J) from below, we use the following elementary estimate. 
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Lemma 5.1. Let p>2 be an integer, and let xi < ■ ■ ■ < Xp be real num- 
bers. Then 



max 

0<r<p \ p — r . r . 

\ J=r+1 J=l 



Proof. For p = 2, the left-hand side just equals For p > 2, 

consider the average x := X]^=2 ^j- Assume first that x < Then 



/ i / *^ 7 I — ^~ / i 

0<r<p\ n — r ^-^ nc ^-^ ^ r) — \ ■> 

\ j=r+l j=l / j=l 



max 



-^{{p-2)x + xi) 
p — I 

p — 2 xi 



2(p-l)^^' + ^^^ p-l 



P / , Xp X\ 
'x„ - Xi) > 



2{p-l 



The case x > ^ii^s boils down to the considered one by the substitution of 
Xj by —Xj. □ 

The definition of b{J) and Lemma 5.1 yield for any block J = + p], 
Plugging this into (5.12) yields that for any k and t, 

,5.13) „„„>.K 0;y{(^)"(^^)"^' 



p=k+l i 



Since the random variables Xj+p — Xj+i , < i < n — p, are identically dis- 
tributed, we get 



p=k+l 
n 

■ n 

p=k+l 



The random variable n{Xp — Xi) is distributed as the sum of p — 1 in- 
dependent random variables having the standard exponential distribution. 



16 



M. LIFSHITS AND Z. SHI 



Therefore, ChernofF's large deviation principle yields that for any fixed 5 > 
and all sufficiently large k and all p> k, 

X„-Xi< = f(x„-Xi< ^{t/T* 

n J \ n 

< pf ""^^^ ~ < (1 + 5){t/T*f 

V P- 1 

< exp{-(l - <5)/((l + <5)(t/r*)2)(p- 1)}, 

from which it follows that 

oo 

nLn{t)>k)<n exp{-(l-5)/((l + 5)(t/r*)2)(p-l)} 

p=k+l 

_ raexp{-(l - 5)I{{1 + 5){t/T*f)k} 
~ 1 - exp{-(l - 5)I{{1 + 5){t/T*Y)} ■ 

We choose now 

i,-u( \ + logra 

''-"^''^^ {l-6)I{{l + 6){t/T*Yy 

so that P(L„(t) > A:) ^ 0, oo. Therefore, 

limsup:^^ < l±4l((l + (^)(t/r*)2)-i in probability. 

n-*oo log n 1 — 

Finally, by sending 5 to zero, we obtain the desired estimate 
lim sup — — < I{{t/T*f r^ in probability. 

n^oo logn 

This completes the proof of Theorem 3.3 for the Poisson model. In order 
to prove it for the i.i.d. -model, we use the representation from Fact 2.1 
with a Poisson process of intensity n. The system with initial positions 
(Xj)j<„ is similar to the Poisson model with initial positions (Xj)j<„, if 
we take similarity coefficients as follows: Cx = {Xn+i)~^, c^, = 
ct = We infer from the similarity that 

(5.14) Lnit) = LnitVx;^), 

where Ln and L„ denote the size of the maximal cluster in the i.i.d. -model 
and in the Poisson model, respectively. Taking into account that Xn+i — > 1 
in probability (law of large numbers), we derive the statement of the theorem 
for the i.i.d. -model from what we have just proved for the Poisson model. 
□ 
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5.2. Proof of Theorem 3.4. Consider first the Poisson model. Let := 
j^*" . The assumptions in the theorem take the form 

(5.15) Hm £n = 0, 

n— »oo 

(5.16) lim ne'i = oo. 

We fix a small f3 £ (0,2) and take a sequence of positive integers (A;„) such 
that 

log(n£^) 
(2-/3)4' 

It follows from (5.13) that 

{L„(t„,)>Mc U [jL{Xi+p-X,+,)<^\. 

It is easier to interpret the condition within the brackets on the right-hand 
side in terms of centered random variables, namely, 

{nXi+p - (i + p)) - {nXi+i + 1)) < -p + l = -pen{2 - e„) + 1. 

Recall that for any n > 1, the random variables {nXi — i)i>o can be viewed 
as sums of independent centered random variables having the standard cen- 
tered exponential distribution. Using the classical Komlos-Major-Tusnady 
estimate [8], we approximate them by a Wiener process W. Namely, we need 
the following: 

Fact 5.2 (KMT construction). It is possible to construct W and (Xi) 
in a common probability space such that for some numerical constants Ci 
and C2, and for all real r > 0, 

(5.17) p( max | {nX, - i) - W{i)\ > rl < (1 + Cin^'^) exp{-C2r}. 

\^l<i<n J 

Since we study only convergence in probability, it is unimportant for us 
whether the probability space in KMT construction is the same for all n or 
not. We just make use of 

f losr TL 1 

lim P<^ max \(nXi - i) - Wii)\ > -|- \ = 0. 

Observe that 

P{L„(t„) > K} 
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n n—p 



f 2 

logn 



Lp=fe„+1 i=0 

+ p( max - - > ^ 

|^l<«<n C2 

The second term on the right-hand side converges to zero by what we have 
just seen. On the other hand, under assumptions (5.15)-(5.16), we have 

P^n > kn£n ~ (2-/3)^^^ ^ logn. Therefore, for any fixed h> 0, 
limsupP{L„(t„) > kn} 

n—foo 

{n n~p ~j 
U \J{\Wii+p)-W{{)\>p£n{2-h)}\. 

By scaling, the probabiHty expression on the right-hand side is 



u u 

'.p=A:„+l i=0 



kn 



"'If + f 



w 



I 

kn 



>^£n{2-h) 

<p( sup sup u-^/'^\W{t + u) -W{t)\>y/Ken{'2-h)\. 



l<u<n/k 

We write, for any a > and 6 > 1 



(5.18) 
Then 



A(a,6):= sup sup u-^/^\W{t + u) - W{t)\. 

0<t<al<u<b 



(5.19)hmsupP{L„(t„) > kn} < hmsupP A — , — > ^/Ken{2 - h) 

n—*oo n— +00 L \ Kn Kn / 

The following lemma gives an upper bound for the tail of A. 



Lemma 5.3. Let A(-,-) be as in [5.18). For any h£ (0,1), there exists 
c/i > such that for all T >1, U >2 and r > 1, 

P(A(r, U)>r)< ChT{log U) exp{-(l - /i)rV2}. 

Proof. By scaling, for any a > 1, 

\W{s + u)-W{s) 



sup sup 

0<s<T a<u<2a 



U 



1/2 



> r =P(A(r/a,2) > r) 



<P(A(r,2) >r). 

Therefore, by the stationarity of the increments of the Wiener process, 

"logC/" 



(5.20) 



P(A(T,C/)>r)< [T] 



log 2 



P(A(1,2) > r). 
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On the other hand, for any centered Gaussian process {Y{v),v G V}, the 
fonowing estimate of large deviations holds (see [9], Chapter 12): 

lim r~^logPl sup > t ) = 

where o"^ := sup^gv'E[y(z;)^]. We apply this estimate to the two-parameter 
process 

r(,,„):=iQi±4^, (,,„)eF:=[0.1]x[l,21, 

SO that lim^_oo?-~2logP(A(l,2) > r) = -i. Therefore, 

P(A(1,2) >r) 
,.>i exp{-(l - /i)r^/2} 

This, in light of (5.20), yields the lemma. □ 
Let now h be so small that 

_ {l-h){2-h)\ ^ 
~ 2(2-/3) 

Applying Lemma 5.3 to T = U = n/kn and r = \fk^{2 — h)en, we get 

< e, "^°g["/^") exp{-(l -h){2- hf elK/2} 

= c,^(n62)-(i+°(i))Mog(n/A:0 
~c,(2-/3)(n4)l-(l+o(l))^ 

(We have used the fact that ne,^ oo.) Since z > 1, we have {ne'^Y~^^~^°^^^'^^ — 
0. Plugging this into (5.19) gives lim^^oo P{-^^n(in) > kn} = 0, which means 
that 



1 > lim sup ■ 

n— >oo fcyi 



= (2 - /3)limsup 

n^oo log(ne^) 

2-/3 L^(tO(r*-t02 . 

= 7^7^ ^™ '^^P 1 — rTF* — r\2V probability. 

{T*Y n^oo \og{n{T* - tnY) 

The desired upper bound for the Poisson model follows by letting /3 — > 0. 
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We proceed now to the proof of the lower bound for the Poisson model. 
We have already established in inequalities (5.6), (5.8) and (5.11) that for 
ah t > 0, 6>0 and all integer k (such that ^ > j^), 

(5.21) F{Ln(,t) < k) <exp{-uPi{k,6)P2{k,5,t)}, 

where u := [j\ and 

Pi(k,6):=w( max f V X,- - - V X,-^ < - + 

P2(A:,<5,t) :=pf^^^ < {1 + 36)-\t/T*f 
\ k + l 

Let us now specify the choice of the parameters. Let /3 > and h > he 
some small numbers. We still denote := rfj" , but this time we set (A;„) 
such that 

log(reg^) 
"'^(2 + /3)4- 

Let Sn ■= hsn- Then 
For all large n, 



{{l + 36n)-\tjT*)^-l)y^k~+T 



> ((1 - 3<5n)(l - En? - l)Vkn + l 



> -i35n + 2en)Vkn + l 
= -{3h + 2)enVkn + l. 
Therefore, for all large n, 

P2(.kn,6n,tn) > p !!^^Ii±l=^^±ll < _ (3/^ + 2)e„ ^^+1 
V V«^n + 1 

Now we use the following result of moderate deviations (see [12], Chapter 
8, page 218). 

Fact 5.4 (Cramer's theorem). Let 1^ be a sequence of i.i.d. centered 
random variables such that Eexp{7|l^|} < co for some 7 > 0. Then for any 
positive sequence (xn) with x„ = o{^/n), we have 
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Since e„ — > [see (5.15)], we are entitled to apply Fact 5.4 to x„ = {3h + 
2)eny/kn + 1. Note that EnVkn + 1 ^ oo [see (5.16)]. Thus, 

P2ikn,6n, tn) > (1 + o(l)) exp{-(3/i + 2)24A:„(1 + o(l))/2} 

(^^ + '^'log(n4)(i + o(l)) 



where z 



(3fe+2)^ 



2{2+/3) 

this choice, (5.16) implies 



> (l + o(l))exp, 
- ^ ^ " 2(2 + /3) 

= (l + o(l))(n4)-(i+o(i))^ 

The parameter h can be chosen so small that z < 1. With 



n 



(5.22) 

Assume for a while that 
(5.23) 



{2 + l3){nei) 
~ 21og(n4) 



2^1-(l+o(l))z 



■ OO. 



liminf Pi(A:„, (5„) > 0. 



Then by (5.21) and (5.22), lim „_^oo P{L„(t„) < kn\ — 0, which implies the 
desired lower bound: 

lim inf > 1 in probability. 



To establish (5.23), we observe that 



(5.24) 



(5„ = h£n ^h(^ 



2 + /3 



2NX 1/2 



OO. 



We use the estimate (5.10), where Qki') is the empirical quantile func- 
tion defined in (5.9). The well-known functional limit theorem for quan- 
tile processes (see [15]) asserts that for A; — > oo, the sequence of processes 

Yk{r) = \fk (Qkir) — r) , r [0, 1], converges weakly to a Brownian bridge W 
in the Skorokhod topology. In light of (5.24), we arrive at 

liminf Pi(/c„,5„) 



> lim inf P I 



max 

n^oo VO<t<l 



> lim inf P( max 

n^oo VO<i<l 



Y^j'iQkAs) - s)ds -^J\q,^{s) - s)ds 



<5r 



1 



max 
Vo<i<i 



l-t 



l-t 

1 o 



1 1 

Xfc„ (s) ds 



tJo 



< 1 



W{s)ds-- fw{s)ds 
t Jo 



< 1 >0. 



Therefore, relation (5.23) is true and Theorem 3.4 is proved for the Poisson 
model. 
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Consider now the i.i.d. -model. We use (5.14) which reduces the problem 
to the Poisson case. Indeed, by the central limit theorem, for any e > 0, there 
exists M > such that 



limsuppf|VX„+i-l| > 

n— >oo V V^/ 



< e. 



Let in ■=tn{l + ^)- The sequence (in) satisfies the assumptions in Theorem 
3.4, and the norming sequences are equivalent: 

j.^ {T*-tn)-^login{T*-tnf) ^ ^ 

(T* - 4)-2 iog(n(r* - 4)2) 

By (5.14), if V^;77< 1 + ^, then L„(4) > Ln{tnVx;^) = Ln{tn). There- 
fore, for any /? > 0, we have, for all large n, 

^~ (T-)^log(n(r*-t„)^) 

^n{tn)> ^2 + /3)(r*-t„)2 



<^[VX;^l>l + —]+¥[Ln{in)> 



M\ , ^fr ,1 ^. (T*)^log(n(r*-y^) 



^) V ^ ' (2 + 2/3)(r*-t„) 

The last probability term on the right-hand side tends to zero by what we 
have proved for the Poisson case. Thus, 

(r*)2iog(n(r*-t„)2) 



limsupP L„(i„) > 



(2 + /j)(r*-t„)2 



M 



< limsupP VXn+i > 1 + ^] <e. 

n—>oo \ V^/ 



Since e is arbitrary, we get the desired upper bound: for any /? > 0, 

The lower bound follows exactly in the same way. □ 

5.3. Proof of Theorem 3.5. We first focus on the lattice model. In this 
model, for every block J = + k], we have 

(5.25) MO) = i E x,(0) = i E ^' + ^ + ' 



Let 



k k n 2n 
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which denotes the average of normahzed initial speeds. Then (4.1) takes the 
form 

2i -\- k -\- 1 

(5.26) x*j{s) = V(TnUjs + -^p{n-2i-k) — , 

Zn 2n 

and (4.5) becomes, for any t>0, 

(/3-a)(l-7pt2) 



(5-27) T*j^f3 < t ^ U(„j] - ^] > 



2nant 



We start with the proof of the lower bound for Ln{t). Let e G (0, 1/3) and 
let k be an arbitrary positive integer. Consider = [{i + e)k, {£+1 — e)k]riZ. 
According to (4.7), for any £, if there exist a < {i + e)k and P > {£ + 1 — e)k 
such that maxjg t* ^ < t, then all the particles of the block belong to 
a common cluster at time t. We apply this to a = £k and P = {£ + l)k to 
see, in light of (5.27), that the condition 

. - , ^ A;(l - 7/Ot') 

mm(n(,,,,.] -n(,. (,+,),]) >^—^ 

ensures a single cluster at time t from the block Ji. Thus, similarly to (5.6), 
we have 

(5.28) P(L„(t) < (1 - 2e)k - 1) < exp{-uPi}, 

where v := [^J denotes the number of blocks and 

k{l - jpt^) 



Pi :=P min(n(oj] - U(^j^k]) > 



Let us express the average speeds in terms of the random walk 

j 



i=l 



Then n(oj-] - UQ^k] = ^%~^) and, thus 



(5.29) 



k I 2nar,t k\ k 



e< ^ < l-e!>. 



Let us now specify the choice of k by letting 



6 = 6„ :=log( ^ ), 



k = kn^{l- Aef'''ci{t){nanf'^h"\ 
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where ci(t) is defined in (3.8). Note that by the definition of ci(t), ^' ) 
(1 — 4e)(26)^/^. Accordingly, for all large n, 

Pi > p|a:-V2 (^Sj - ^-Sk^ > (1 - 3e)(26)V2^ (^1 _ , ^ < ^ < i _ ^ 
Let W be the Wiener process in (5.17). We have 
Pi > p|fc-V2 (^^(^-^ _ lw{k)^ > (1 - e)(26)V2^ (^1 _ £ 

k-'/'\W{j) - 5,1 < £(26)1/2 1. (^1 _ ^ ^ < J < 

> F{k^^/'^{W{ks) - sW{k)) > (1 - e)(26)i/2^(l -s), e<s<l-e, 



max |Vr(j) - Sj\ < e\l - e){2by/A. 

Since s i-^ ^"-'^/^^(/j;^)^ s > 0, is also a Wiener process, we get Pi > P2 — P3, 
where 

P2 := P{VF(s) - sW{l) > (1 - e)(26)i/2s(l - s), e < s < 1 - e}, 
P3 := p( max |W(j) - 5,| > £^(1 - e){2bk)^/^}. 
In order to estimate P2, we apply the functional large deviation princi- 

o 

pie (see, e.g., [9], Chapter 12, Theorem 6) to the Brownian bridge W{s) = 
W{s)-sW{l): 

liminf 4TlogP{P7(s) - sW(l) >Rs(l-s), e<s<l-e} 

R^oo P^ 

By (3.5), we have b 00, which implies that for large n, 
(5.30) P2>exp{-(l + e)2(l-e)26/3}: 



n \ (i-s')V3 



We now show that P3 is negligible compared to P2. Indeed, by the Komlos- 
Major-Tusnady estimate (5.17), 

(5.31) P3 < (1 + Ci/fci/^)exp{-C2e^(l - e){2hk)^/^}. 

Since | ~ (1 - 4e)2/3ci(t)(^)2/3 > 1^ we have, for all sufficiently large n, 

C2e^(l - e)(26/c)i/2 > ^ + log(l + Cik^^"^), which implies that for all large n, 

-1 



P3 < exp{-6} 



AGGREGATION IN SELF-GRAVITATING GAS 25 
In view of (5.30), this yields P3. Hence, 

uPi > U{P2 - P3) 



k^^ - ci(t)(na„)2/3(log(n/(72))V3 \al 



The expression on the right-hand side equals ^7Ty(^)'"'^ ^ l/^(log(^)) 
and thus goes to infinity. In light of (5.28), we obtain 

liminf ,2/3n"? / 2^M/3 ^ " " 4e)'^'ci(t) in probability. 

The lower bound in Theorem 3.5 for the lattice model is proved, since e can 
be arbitrarily small. 

We proceed now to the proof of the upper bound (for the lattice model) . 
We have already seen in the proof of Theorem 3.3 that if L^it) > k, then 
there exists a block J = [i^i + p] of size p between k + \ and 2k, that is free 

J and such that t j 



up to its collapse time t j and such that fj <t. Let, for simplicity, p be even 



say, p = 2q. By virtue of Proposition 4.4, inequality fj < t means that 



In particular, we have T*j_,_g j_,_2q < t, which, in light of (5.27), is equivalent 
to 

o(l-7pt^ 

Similarly, whenever p is odd, say, p = 2q + 1, we arrive at 



'^(iA+q] - '"(i+g+l,i+2g+l] > 



(1 - ipt") 

na„t 



Summarizing, we have 



{L„(t)>/c}cOiU02, 



where 



k n—2g 



and 

k n—2q 



n 



2 ■ 



II IN ^ \-( ,^ q''\i-ipt^) 

U U i ^ Z^("i+i - Ui+q+i+j) > y= 7 
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Thus, 



k 

F{Ln{t) >k}<2n J2 Pg^ 

q=k/2 



where 



1 " 



and 



Let e > 0, and let 

k = k{n) ~ (1 + 2e)c2(t)(n^7„)2/3(log(na2/5))i/3^ 

where C2(t) is defined in (3.8). Since we want to apply Fact 5.4 to Pq, let us 
check that 

(5.32) l™%| = 

uniformly over q < 2k{n). Indeed, 

gV2 \nanJ \\ nan J 
By using twice the assumption (3.5), we have 

log(ncr^^/^) 21og(n) 
iim sup < lim sup = U 

n—*oo nCn n—*oo nCn 

and (5.32) follows. Now, by Fact 5.4, for all large n and all q G [k/2, k], 
^ ^^n^(TT7) j = ^"n-2(l + .) 2(n.„)^t^ j' 

so that 
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The expression on the right-hand side is, when n is large. 



(log(na^/'))2/3 "^l 25*2 



.(log(na^/'))2/3 
which goes to as n ^ co. Hence, 

hmsup ^^^^ — < (1 + 2e)c2(t) in probabihty. 

(na„)2/3(log(na^/'))i/3 

This yields the desired upper bound in Theorem 3.5 for the lattice model. 

The theorem is thus proved for the lattice model. It turns out that in the 
considered range of initial speeds, the fluctuations in initial positions have 
no significant influence upon the asymptotics of Ln{t). Indeed, for arbitrary 
initial positions, (5.29) becomes 



(5.33) 



where 



^<|<l-£|, 



3 k 
Dj ■■= -iZi^^Xi -i)- V {nX, - i). 

Hence, for any ei > 0, we have Pi > Pi — Pa, where 

Pi :=P(^mm(n(o,] -^o-,,]) > j, 

P4:=pfmax|Z^,|>^iMl^^). 
V 2 / 



In the Poisson model, it follows from the large deviation principle that for 
some constant c = c{'y, p,t,ei) G (0,oo), 

(5.34) P^<exp{-ck}, 

whereas from what we have proved for the lattice model, we know that 
Pi ^> exp{— cfc}. Therefore, the proof of the lower bound for the Poisson 
model goes through along the same lines as for the lattice model, with Pi 
in place of Pi . 
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The same happens with the upper bound, where Pq should be replaced 

by 

Again, the exponential bound in (5.34) suffices to conclude the proof as for 
the lattice model. 

The passage between the Poisson and the i.i.d. models via Fact 2.1 is 
straightforward. 

5.4. Proof of Proposition 3.6 (a sketch). Few changes are needed with 
respect to the proof of Theorem 3.5, except that we have to provide alterna- 
tive tools to those based on exponential moments. For the lower bound, we 
can replace the Komlos-Major-Tusnady estimate (5.31) by the Sakhanenko 
bound (see [13]) which states that under condition of the finiteness of the 
pth moment, 

P max >r < VfcGN, Vr>0, 

where C > is an unimportant constant. For the upper bound, we can use, 
instead of Cramer's moderate deviation principle (Fact 5.4) in the estimate 
of Pq, the following result of Nagaev [11]: for a G (0,1) and i.i.d. random 
variables Yi, . . . ,Yq with mean zero, unit variance and finite pth moment, 

IP^E^i ^ vVQ^ ^ CaAl^'^^^nYl^y" + exp{-ayV2}] 

VgGN, Vy>0, 

where Ca,p > is an unimportant constant. The rest of the proof is along 
the same lines. 

5.5. Proof of Theorem 3.7 (a sketch). Under (3.10), the fluctuations of 
speeds and initial positions have comparable influence on the asymptotic 
behavior of L„(t), hence, we have to take both of them into account. Fortu- 
nately, the large deviation principle provides less precise (than in the preced- 
ing proofs), but still sufficient, estimates. We just outline how the estimation 
works for the upper bound for the Poisson model. The changes in the proof 
of the lower bound are similar. 

Fix M > 0; let g > Mlogn and B := i^^. We obtain for Pq in (5.35) 
that for any e > and all large n. 
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Mq{l -£)(!- -fpt^) 
ct 4 



) 



< g-(l+o(l))/i(B)g ^ g_{l+o(l))/2(C)<7 




q>Mlogn 



and the upper bound follows. 
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